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1. Introduction 

For the era of precision LHC physics(l% or better total theoretcial precision []]]]), we need [||, 
|, § resummed @(o%L n ),ff(a s aL n '), (ff(a 2 L n ") corrections for n = 0, 1,2, ri = 0, 1,2, n" = 1,2, 
in the presence of parton showers, on an event-by-event basis, without double counting and with 
exact phase space. We present the first step in realizing our new MC event generator approach to 
such precision LHC physics with amplitude-based QED0QCD resummation []|] by introducing 
the attendant new parton shower MC for QCD that follows from our approach. This will set the 
stage for the complete implementation of the QED&QCD resummed theory in which all IR singu- 
larities are canceled to all orders in a s and a. This new parton shower MC, which is developed in 
the HERWIG6.5 [g] environment and which we have called HERWIRI 1.0(31) [0], already shows 
improvement in comparison with the FNAL soft pr data on single Z production as we quantify 
below. We also note that, while the explicit IR cut-offs in the HERWIG6.5 environment will not be 
removed here, HERWIRI only involves integrable distributions so that in principle these cut-offs 
could be removed. 



We first review our approach to resummation and its relationship to those in Refs. Sec- 
tion 3 contains a summary of the attendant new IR-improved DGLAP-CS [1C, 11] theory [|l2|]. 
Section 4 presents the implementation of the new IR-improved kernels in the framework of HER- 
WIG6.5 to arrive at the new, IR-improved parton shower MC HERWIRI 1.0(31). We illustrate 
the effects of the IR-improvement with the specific single Z production process at LHC energies. 
We compare with recent data from FNAL to make direct contact with observation 1 . 



2. Review of QED^QCD Resummation 

In Refs. [jl, [l2]] we have derived the following expression for the hard cross sections in the SM 
SU 2L x Ui x SU§ EW-QCD theory 

,p. _ sum ir (QCED) v 1 / ^Pi d 3 q2 A- d i k jl jL d 3 k' j2 
£XP J^n\m\J q o 11 k . 11 Vj% 

_gy e tr{px+*-Pi-*-T* h -U' J2 )+d Q ce» ~p (jt . ^ . ^ (2 .i) 
[2nf 

where the new YFS -style [ ]14[ ] residuals j8 n m (^i, . . . ,k n ',k^, . . . ,k' m ) have n hard gluons and m hard 
photons and we illustrate the generic 2f final state with momenta p 2 , <?2 specified for definiteness. 



The infrared functions SUMir(QCED), Dqced are defined in Refs. [Q, |12p. Eq. (2.1) is exact to 
all orders in a and in a s . 



Given that the approaches in Refs. M, W to QCD resummation have been shown in Refs. [15] 



to be are equivalent and given that we show in Refs. [ ]12p that our approach is equivalent to that in 
Ref. [§], it follows that our approach is also equivalent to that in Ref. See Refs. [|l2|] for the 
attendant further details. 



'From Ref. | Jl3| ] the current state-of-the-art theoretical precision tag on single Z production at the LHC is (4.1 ± 
0.3)% = (1.51 ±0.75)%(QCD) ®3.79(PDF) 0O.38±O.26(£W)% and the analogous estimate for single W production 
is ~ 5.7%. 
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3. Review of IR-Improved DGLAP-CS Theory 



The result Eq. ( 2. 1 ) allows us to improve [ 12] in the IR regime the kernels in DGLAP-CS [ 10, 
11] theory as follows, using a standard notation: 



PZ P ^=C F F YFS ( Yq )e^ 



l-z 



{l-zy"-f q {y q )8{\-z) 



PG*(z)=C F F YFS ( 7ll )e^ 1 + {1 Z)2 



P^(z)=2C G F YFS ( 7G )e^{ 



l-z 



7c 



+ t^(1 -zV a 
1 — z 



+ 2^ 1+rG (l -z) + z(l -z) l+ K)-f G (r G )8(l -z)}, 



YFS 



(YG)e± Sa Uz 2 (l-z) ra + (l~z) 2 zK}, 



(3.1) 



where the superscript "exp" indicates that the kernel has been resummed as predicted by Eq. (2.1) 



when it is restricted to QCD alone and where we refer the reader to Refs. Ql 2fJ for the detailed 
definitions of the respective resummation functions Fyfs, 7a, 8a,/g,A = q,G 2 . 



We refer the reader to Refs. [12, 18] for discussion of a number of illustrative results and 



implications of the new kernels and Eq. (2.1) that are beyond the scope we have here. The net 
effect of the results in Refs. [|l2|, 18] is that we have a consistent theoretical paradigm based on 



Eqs. ( |2. ljp. 1| ) for precision LHC theory that can be systematically improved order-by-order in 
perturbation theory with no double counting. With an eye toward the full MC implementation of 
our approach, we turn next to the initial stage of that implementation - that of the new kernels. 



4. Realization of IR-Improved DGLAP-CS Theory via MC Methods 



We have implemented the new IR-improved kernels in the HERWIG6.5 environment to pro- 
duce a new MC, HERWIRI1.0, which stands for "high energy radiation with IR improvement"] 19]. 



Specifically, we modify the kernels in the HERWIG6.5 module HWBRAN and in the atten 
dant related modules 3 with the following substitutions: DGLAP-CS P AB => IR-I DGLAP-CS P exp 



AB 



while leaving the hard processes alone for the moment. We have in progress [20] the inclusion in 



our framework of YFS synthesized electro weak modules from Refs. [21]for HERWIG6.5, HER- 
WIG++ [0], and MC@NLO ^ hard processes 4 , as the CTEQ [|2§] and MRST(MSTW) [^]] best 
(after 2007) parton densities do not include the precision electroweak higher order corrections that 
are needed in a 1 % precison tag budget for processes such as single heavy gauge boson production 
in the LHC environment 

The details of the implementation are given in Refs. [f7|, [H]] and we do not reproduce them here 
due to a lack of space. We have done many comparisons of the properties of the parton showers 
from HERWIG6.510 and HERWIRI1.031. In general, as we show here in Fig. [j], the IR-improved 
showers tend to be softer in the energy fraction variable z = E/Eseam where E (Eseam) is the cms 



2 The improvement in Eq. (3.1) should be distinguished from the also-important resummation in parton density evo- 
lution for the "z —t 0" regime, where Regge asymptotics obtain - see for example Ref. [|l6[ |l7|. This latter improvement 
must also be taken into account for precision LHC predictions. 

3 We thank M. Seymour and B. Webber for helpful discission. 

4 Similar results for PYTHIA [ p4| and for the new kernel evolution in Ref. [ p5| ] are under study. 
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Energy fraction distribution of parton shower for single Z production. 
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Figure 1: The z-distribution(ISR parton energy fraction) shower comparison in HERWIG6.5. 



parton(beam) energy for hadron-hadron scattering respectively. See Refs. [|7], [TJ] for the complete 
discussion of such comparisons. We have also made comparison analyses with the data from FNAL 



on the Z rapidity and pj spectra as reported in Refs. [|28J, [290. We show these results, for 1.96TeV 
cms energy, m Fig. g. We see that HERWIRI 1.0(31) and HERWIG6.5 both give a reasonable over- 
all representation of the CDF rapidity data but that HERWIRI 1.031 is somewhat closer to the data 
for small values of Y. The two £ 2 /d.o.f are 1.77 and 1.54 for HERWIG6.5 and HERWIRI 1.0(31) 
respectively. The data errors in Fig. |2|(a) do not include luminosity and PDF errors [28], so that 
they can only be used conditionally at this point. We note as well that including the NLO contri- 
butions to the hard process via MC@NLO/HERWIG6.510 and MC@NLO/HERWIRI1.031[j23|] 5 
improves the agreement for both HERWIG6.510 and for HERWIRI1.031, where the # 2 /d.o.f are 
changed to 1.40 and 1.42 respectively. That they are both consistent with one another and within 
10% of the data in the low Y region is fully consistent with what we expect given our comments 
about the errors and the generic accuracy of an NLO correction in QCD. A more precise discus- 
sion at the NNLO level with DGLAP-CS IR-improvement and a more complete discussion of the 
errors will appear [30]. These rapidity comparisons are then important cross-checks on our work. 
We also see that HERWIRI1.031 gives a better fit to the DO p T data compared to HERWIG6.510 
for low pj, (for pj < 12.5GeV, the # 2 /d.o.f. are ~ 2.5 and 3.3 respectively if we add the sta- 
tistical and systematic errors), showing that the IR-improvement makes a better representation of 
QCD in the soft regime for a given fixed order in perturbation theory. We have also added the 
results of MC@NLO [23] for the two programs and we see that the &(cc s ) correction improves 
the # 2 /d.o.f for the HERWIRI1.031 in both the soft and hard regimes and it improves the HER- 
WIG6.510 # 2 /d.o.f for pj near 3.75 GeV where the distribution peaks. For pj < 7.5GeV the 
£ 2 /d.o.f for the MC@NLO/HERWIRI 1.031 is 1.5 whereas that for MC@NLO/HERWIG6.510 is 
somewhat worse. These results are of course still subject to tuning as we indicated above. We await 
further tests of the new approach, both at FNAL and at LHC. 



We thank S. Frixione for helpful discussions with this implementation. 
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Figure 2: Comparison with FNAL data: (a), CDF rapidity data on (Z/y*) production to e + e pairs, the 
circular dots are the data, the green(blue) lines are HERWIG6.510(HERWIRI1.031); (b), DO p T spectrum 
data on (Z/y*) production to e + e~ pairs, the circular dots are the data, the blue triangles are HERWIRI1.031, 
the green triangles are HERWIG6.510- in both (a) and (b) the blue squares are MC@NLO/HERWIRI 1.031, 
and the green squares are MC@NLO/HERWIG6.510. These are untuned theoretical results. 
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